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We develop a method for solving the field equations of a quadratic gravitational theory coupled to 
matter. The quadratic terms are written as a function of the matter stress tensor and its derivatives 
in such a way to have, order by order, a set of Einstein field equations with an effective T^v We 
study the cosmological scenario recovering the de Sitter exact solution, and the first order (in the 
coupling constants a and (3 appearing in the gravitational Lagrangian) solution to the gauge cosmic 
string metric and the charged black hole. For this last solution we discuss the consequences on the 
thermodynamics of black holes, and in particular, the entropy - area relation which gets additional 
terms to the usual jA value. 
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a^ 

pH . I. INTRODUCTION 

^ . 

Higher derivatives gravitational theories have been proposed at classical level as an extension of Einstein's theory 
in an attempt to unify other fields with gravity |lj and to avoid the cosmological singularity Ol . Quadratic curvature 
counterterms also appear in the renormahzation of the one- loop semiclassical approximation p| . The pure Quadratic 
theory is renormalizable Q and asymptotically free H although it has problems with unitarity Q. Higher order 
^ . gravitational theories can lead naturally to inflation HR] and arise as the low energy limit of string theory p] . 
^/^ ' For the sake of definiteness we deal with the following Lagrangian formulation of quadratic theories 

S = SG + Sm^ I d'^xy^l - 2A + i? + ai?2 

o. 

^ . +l3Rt..R'"' + 2A:/:,n \ , (1) 

G^ 

^ i' where we have dropped the R^^xp term by use of the Gauss-Bonnet invariant in four dimensions |^ . 

I . The field equations derived by extremizing the action S are given by (we use the sign conventions of Rcf. pX|) 

5-1 . 

bJO. 1 

Rfiiy - -^RQiiu + ^9iiu + aH^^ + I31^j,y 

= kT,.^-^^, (2) 



in 



X 



where we have chosen units such that c = 1 and k = IGttG, and where 

1 



and 



2 
Note that the trace of equation (0) takes the simple form 



H/j^i, — —2R-^i, + 2g^,yBR — —gp,i,R + 2i?i?^,^ , (3) 



Itiw - ~^R(^ ■„)a + ^Rf^iy + l^dtJ^i^^R 

+2R"Rai, — -gfii^RapR" ■ (4) 
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2(3a + /3)Di?-i? + 4A = fcT . (5) 

In section 2 of this paper we develop a method to find a metric solution of the classical field equations of higher 
order gravitational theories with sources by successive perturbations around a solution to Einstein Gravity, which 
will represent for us the zeroth order. In section 3 we apply this method to find first order solutions in the coupling 
constants a and /3 for the straight gauge cosmic string and the charged black hole. For this latter solution the 
associated thermodynamics is studied and the corrections to the Bekenstein - Hawking temperature and entropy are 
discussed. We end the paper with some further discussion on this first order solutions. 

II. PERTURBATIVE SOLUTION 

From the Eq. (||) we can see that R satisfies a massive scalar wave equation. For this field to have a real mass we 
impose 

3a + /3 > , /? < (6) 

where the last inequality comes from linearizing Eq. ^ and asking a real mass for the massive spin-two field i/'/ji' 
related to R^^j, (see Ref. 0). Conditions (^ are called the no-tachyon constraints. 

The coupling constants a and /3 have to be at most of the atomic scale, since otherwise they could have observable 
effects in, for instance, the solar system or binary pulsars. In the Ref. Il2] it was found that a < 10^^^ by requiring 
that the inflationary period be of sufficient duration (where Ipi is of the order of 10"'^^ cm). A similar bound for f3 
can thus be obtained. 

We will then consider only small curvatures in our method, such that 

a|i?|<l , |/3i?^,|<l. (7) 

Thus, we can obtain perturbative solutions to the field equations in a power series of a and f3. We will take as 
the starting metric, gliJ , a solution of Einstein equations and then, systematically, successive orders of the coupling 
constant a and p. 
To zeroth order 

R,A0 - iR^'^g'-^} + A<?(°J = fcT,.(ffW) , (8) 

with trace 

i?(o) = 4A - fcr(o) (9) 

Once we solve for the Einstein metric gfiJ , we pursue the iteration to the first order 

i?(i)_ii^(i)„(i)+AoW 

= fc^M^ {g^} ) - o^H^, (g(,°J ) - (31^. {gf} ) , (10) 

where to first order_ in a and /3 it is enough to consider in Eq. (|To|) Hf^J and I(i,J . 



By use of Eqs. (g) the right hand side of Eq. ( |10[) can be written as 

T(l)e// 






where 

!<)- ^i/..(3(°\T(«)) =. 2r(;) - 25(°jnT(°) 



k 
k 
2'- 



+ ^ffl°.HT(0))2 _ 2fcr(0)r(0) + SATW - 2AT(°)g,. , (12) 



and 



Thus Tf^iu plays the role of an effective energy momentum tensor in an Einsteinian field equation for g/ji/. It is 

easy to see that T°^. '^ — Q, i.e. satisfies a conservation law (with respect to .g),J). T°^ also inherits the symmetry 
properties of T^^. This properties will hold to every order of the approximation. 
We can now generalize to the n-th order approximation 






(14) 



Expressions (O)-(O) greatly simplify when T^^ is diagonal and its components depend essentially on only one 
coordinate; let us call it r. This will be the case in the applications we will deal with in the next section (we will also 
take A = for simplicity). Thus, in this case 






(15) 



and 






gfifi 
+g"X 



r ,rr ' ra \ r ,r a ,r ) ) ^ fi MM M i 

2A;Tr; - s'''^ v.- + .9""r:, jr; ,, - t,) 



r. - ^T^ - 2fc(T/)2 + ^T'^'^T^^ 



(16) 



where the metric dependence in this expressions is with respect to the zeroth order, i.e. solution of usual Einstein 
equations, Eq. (||). In Eqs. ( [15|) and (|lq ) sum is over a index but not over fj,. 



III. FIRST ORDER SOLUTIONS 

We are now ready to compute the different metrics solution of the fourth order field equations. Three main 
astrophysical scenarios where gravity plays an important role can be studied: cosmology, topological defects and 
black holes. 



A. De Sitter Universe 



It is a solution of vacuum Einstein equations with cosmological constant. This solution is useful for describing 
an inflationary phase in the very early universe where higher order gravity terms are presumably worth taking into 
account. 

In our approximation, for T^j, = 0, we have that de Sitter metric is a solution of the fourth order field equations order 
by order (see Eq. (Oh), without restrictions for the value of A. Thus, it is an exact solution as can be directly verified 



from the field equations |1^. In fact, in general, vacuum solutions to Einstein equations (even with cosmological 
constant), are solutions to the quadratic theory (the converse, in general, is not true). 

The Robertson - Walker metric can be also replaced in Eqs. (|l^) and (|l|) to find the corrections to the general 
relativistic results. However, this approximation breaks down near the cosmological singularity and we already dispose 
of some exact solutions (see Ref . H ) . 



B. Gauge Cosmic Strings 



Our approximation is specially suitable for dealing with the gravitational field of topological defects since while its 
effects are of relevance they are relatively small compared to to the planck scale, thus a first order approximation in 
a and /3 should provide acceptable results. 

Topological defects are expected to be formed during phase transitions driven by the evolving (and cooling down) 
early universe whenever the manifold of equivalent vacua after the spontaneous symmetry breaking is not shrinkable 
to a point. We can model this process by studying an n-internal components scalar field with a Mexican-hat like 
effective potential with an associated energy scale of symmetry breaking, 77^, of order 10~^ for the GUT scale Il3] 
and coupling constant A. Depending on the dimension of the vacua manifold the topological defects can be domain 
walls, cosmic strings or monopoles and depending on whether the (internal) symmetry that breaks down is a local 
(or gauge) or a global one, the formed topological defects will be localized in a tiny core of size Tc ~ —K^ or infinitely 

extended 0. 

To see explicitly the effects of this higher derivative theory of gravity we chose one particular topological defect. 
Since in the next section we will deal with a spherically symmetric system (a charged black hole) we will deal here 
with the local (or gauge) straight cosmic string that posses cylindrical symmetry. The other topological defects can 
be treated in an analogous way (see Ref. lO] for a thorough account of them). 

A straight, static, cylindrically symmetric local string lying along the z — axis can be characterized by the following 
energy-momentum tensor (for r ^ re) 



V 



T, 



Tf'5{r) 



2TTry/B{r) 



T 



-T, 







For a generic metric of the form 



ds^ ^ A{r){-dt^ + dz^) + dr^ + r^B{r)d0^ , 
the exact metric, solution of Einstein equations is given by p6| . 

?7^ 



Air) = 1 , B{r) = Bq - U 



(17) 



(18) 



(19) 



By plugging the generic metric (lia) into Eqs. (nfl) and (O) we can now write the solution to our first order Einstein 
equation with effective source T°^ 



A{r) ^ci + k 



r r"dr" r ^[Tf<'-Tf-] 



(20) 



and 



B{r) = C2 - ^^ / dr" 



dr r 



-T-{3a + l3)DT +A{r") 




TS 



(21) 



where ci and C2 are arbitrary constants to be conveniently chosen in such a way to recover general relativistic results. 
Explicitly replacing here the form of the gauge cosmic string stress tensor given by Eq. ( |l7|) and upon integration 
and regularization of delta squared terms, we obtain the metric components up to linear terms in a and /3 



A{r) = l + {2a + (3)kT 



(22) 



2 \ 2 ^ 2 \ ^ 

B(r)=\\-k^\ -4a/c(l-fcyj T 



(2^)4 



(23) 



Wc observe here that the corrections to the general relativistic metric due to including quadratic terms in the curvature 
in the gravitational lagrangian can be classified in two types. The terms in v4(r) and -B(r) proportional to the trace 
of the matter stress tensor, T, give localized contributions. They are different from zero only in the core of the gauge 
string and vanish for r > r^- This is essentially what was found in Ref. |ll[| , where we restricted the analysis to 
linearized terms in the curvature and only considered its higher derivatives. It is precisely the additional terms not 
considered in Ref. |11| , i. e. quadratic in the stress tensor, that generate the non-localized term, proportional to r^"^ 
appearing in B{r). The structure of this term is such that it is linear in the coupling constants a and (3 (due to our 
approximation), and if one considers they have an associated radius ri, the dependence (ri/r)^ is the only extended 
possible one not divergent as r ^ cxd. On the other hand the factor k^r]'^ already appears in processes such as particle 
production by the formation of Cosmic String ||l^ and Global Monopoles [Q . We see that due to the quadratic gravity 
terms, space outside a straight local cosmic string is no longer flat as in General Relativity, but curved with curvature 
terms typically going as r~^. This dependence also appears when one considers the renormalized energy-momentum 
tensor due to vacuum polarization po|-p2|. 

C. Charged Black Holes 

We shall now study spherically symmetric solutions to the quadratic field equations in the first order approximation 
Eqs. (|lO|)-(|ll|), which represent charged black holes; starting from the General Relativistic solution, i.e. the Reissner- 
Nordstrom metric, 



ds2 ^ gtt{r)dt^ + grr{r)dr^ + r^dV? , (24) 

where -gtt{r) = 5rr(r)"^ = (1 - 2M/r + Q'^jr^) and d^^ == d-d'^ + sin^ ^dip^ . 
The non- vanishing components of the energy-momentum tensor are |1G| 

Tl = T; = -Tt = -T; = -^ . (25) 

The exact metric, solution of the Einstein equations with source, can be written as |23| ] (in the Schwarzschild gauge, 
Eq. (H)) 



g-}(r) = \ + - / f'Tf'df 



(26) 



attir) = -grrirY^ exp U^ (T; - Tlf^fg„{f)df\ . (27) 

By use of Eq. ( psf ) to compute the T^JJ (Eq. (|ll|)) we find the first order corrections to the energy-momentum tensor 

_4/3Q2 2Af Q2 

-^;:^(i- — + 7^)' (28) 

and replacing it into Eqs. (^) - ( p7| ) we obtain 

1 2Af g2 12/3Q2 (\ M Q'\ 

- gttir) ~ grrir) e~^^ 



2M 
r 



Q2 2/?Q2 



1- 



M , Q2 



(30) 



We observe here that the a couphng constant does not appear. This is due to the fact that the trace of the electro- 
magnetic energy momentum is zero. In fact, Whitt |Q has shown that for the quadratic theories coupled only to the 
a-term (i.e /3 = in Eq. (g)) there exist a "no hair" theorem stating that the only black hole solution (with spherical 
symmetry) , must be the Reissner - Nordstrom family. This can be directly seen from our method, since order- by-order 
the a-contributions vanish, thus leaving us with the Reissner - Nordstrom solution. This is not the case, of course, 
when one allows the (3 terms be different from zero, and expressions Eq. ( pO[ ) represent the 0{P) generalized Reissner 
- Nordstrom metric. 

Direct inspection of metric (pfl) allow us to see the improved result that generates the perturbative method presented 
in this paper with respect to the "linearized" approach of Ref. [p5| . There, only terms up to 0(r~*) have been 
considered. 

The validity of this metric will be assured if the condition Eq. (0) holds. In our case, this takes the form —fi/r^ <C 1. 
Thus, Eq. (pO|) will be a good approximation to the charged black hole solutions in quadratic theories if 



TH > 



-P 



(31) 



where th is the radial coordinate of the event horizon (in Schwarzschild's gauge). 

The geodesic motion acquires an additional small repulsive term in the effective potential since for the no-tachyon 
constraint /3 < 0. This provides, in principle, a way of detecting the physical effects produced by the higher order 
corrections to the gravitational lagrangian. 

The radial coordinate of the horizon can be computed directly from making vanish gttifH) given by Eq. (p3C 



TH 



02 '1 



3Q^ 



1 - 



r+=M + V^ 



(32) 



We observe that the quadratic black hole shrinks with respect to the corresponding general relativistic one. 

The extreme black hole will be now reached with a maximal charge lower than the general relativistic one (/? < 0) 



Q 



2 

max 



5 



(33) 



thus th given by Eq. (p2), will remain always bounded. In fact, r^* 
The horizon area, Aj/, will then be 



M 



JL 

5J\/ 



A 



H 



inri 



2pQ'' 



i'-m 



1- 



(34) 



which is smaller than in the general relativistic case. 

The electric potential on the horizon can be independently computed by use of Maxwell equations in curved 
spacetime Eq. (|30|) 



^H(r^ 



Qdr 



H r^\/-9ttgr 



Q_ 



2/3Q2 /3-2QVr 



5r1 



1 - Q2/r 



2/^2 



(35) 



Let us now turn to the thermodynamical properties of black holes and see how they are modified by the quadratic 
corrections. 

We can easily compute the Bekenstein-Hawking temperature from the surface gravity, k: 



T, 



1 



H 



9tt 



27r 



47r y/~gttgr 



(36) 



Thus, 



T, 



H 



Airru 
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/?Q4 /2-QVr 



STrrT 



1 - Q2/^4 



(37) 



Since /? < 0, we observe that the effect of the quadratic gravitational theory corrections will be that of decreasing the 
black hole radiation temperature with respect to the general relativistic value (with the same M and Q). Leaving 
thus out open the possibility of switching off black hole evaporation and leaving behind a charged remnant with a 
mass of the order of the Planck mass. 

By inverting Eq. ( p4| ) for M as a function of Ah and Q. We obtain the fundamental relation of black hole 
thermodynamics. Differentiation of this equation produces the first law, 



dM ^ 



dM 



Q dQ 



dQ = TadS + ^„dQ , 



(38) 



that can be used to obtain the entropy of the quadratic black hole, S. 
Since, 



dM 
'dA 



1 
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Q 167rr+ V rl 






and 



dM 
'dQ 



Q oQ 



6Q" 

5 r-f 



(39) 



3g 

5 r7 



(-f) 



(40) 



we can obtain the entropy of the black hole given <^h and Tjj by Eqs. (p3^(p 

87r2/3Q2 



A 



A 



(41) 



We observe here that the entropy of the charged black hole is not simply one quarter of the area, but in general will 
be a more complicated function of A. The effect of the corrections linear in (3 will then be that of increasing the 
gravitational entropy by an amount proportional to Q"^ j A. We can also see that the other parameters characterizing 
the black hole, such as the charge Q (and the angular momentum if we had considered rotation (see Ref. |25|)) enter 
explicitly in the equation defining the entropy. 



IV. DISCUSSION 



In this paper we have presented a perturbative approach to find solutions to the classical field equations of a 
quadratic (in the curvature) theory of gravitation. Since this theory can be considered a generalization of general 
relativity, we start the solutions from the general relativistic metric and then add corrections (presumably small) of 
successive order in the coupling constants a and /3. It is worth to stress here that our method does not obtain all 
the possible solutions to the quadratic theories, but only those expandable around a metric, solution of the Einstein 
equations, in powers of the coupling constants a and /3. We have thus analyzed the three possible scenarios of 
application: Cosmology, Topological Defects and Black Holes, and obtained the first order corrections to the general 
relativistic results. Our method besides allows to make a systematic study of the higher order corrections by means 
of symbolic computing programs |26| . 

The corrected metric (p3) may be of relevance to study the evolution of gauge strings in the early universe for 
both the structure formation scenario where the r~^ term in S(r) could play an important role and for the collision 
simulations where the short range contributions may change the predictions of a string network. A detailed study of 
these effects might provide a link between observation and the a and /3 parameters. 



The charged black hole metric in quadratic theories gets only modified by the (3 coupling. The Reissner-Nordstrom 
metric is no longer a solution to the problem and one must study a different solution (Eq. (|30|)). The thermodynamics 
of this charged quadratic black holes is different from that of a Reissner - Nordstrom black hole. The Bekenstein- 
Hawking temperature as well as the other thermodynamical parameters acquire corrections with respect to its general 
relativistic values (in Ref. [E5[ they coincide with those of general relativity) as one would expect for a different 
theory of gravitation (see Ref. ||2^ and Ref. |g8| for a discussion in the context of string and Kaluza-klein theories 
respectively). 

It is also interesting to remark that for the black hole solution, the relation between entropy and area is no longer 
the simple S = jA, but that given by Eq. {nW . we stress that this modification of the entropy - area relation have been 
obtained by considering a four dimensional quadratic gravitational theory minimally coupled to the electromagnetic 
field. We also expect that in an exact black hole solution this simple equation will break down leaving its place to a 
more fundamental relation p9| , p0[ . 
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